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RAY  THEORY  VS  THE  PARABOLIC  EQUATION  IN 
A  LONG-RANGE  DUCTED  ENVIRONMENT 


INTRODUCTION 

An  important  type  of  underwater  acoustic  environment  involves  a  deep  water  SOFAR  channel 
with  a  superimposed  surface  duct  capable  of  trapping  acoustic  energy.  While  considered  to  be  most 
representative  of  winter  conditions  (the  profiles  contain  no  surface  heating  effects),  this  environment 
prevails  over  northern  latitude  areas  throughout  the  year  and  is  adequate  for  the  intended  parameter 
study.  In  this  report,  we  examine  the  effects  on  acoustic  propagation  predicted  by  a  computer  model  of 
three  ducted  sound-speed  profiles  for  a  nominal  frequency  of  300  Hz.  The  purpose  of  this  study  was  to 
better  understand  the  capabilities  and  limitations  (both  theoretical  and  numerical)  of  a  ray-theoretic 
model  for  predicting  mid-  and  long-range  signal  transmission  loss  (TL)  in  strong  surface  ducts,  for 
example  as  in  the  winter  North  Atlantic.  Several  source  and  receiver  configurations  (combinations  with 
one  or  both  in  and  out  of  the  duct)  have  been  studied.  To  test  validity,  the  ray-theoretic  calculations 
are  compared  to  previously  accepted  results  predicted  by  a  propagation  model  based  upon  the  parabolic- 
equation  (PE).  Both  model  predictions  agree  to  first  order  (TL  measurements  and  acoustic  field  pat¬ 
terns  are  comparable),  while  any  second-order  discrepancies  are  explained  by  the  theoretical  or  numeri¬ 
cal  limitations  of  the  different  approaches.  We  conclude  that  when  both  the  source  and  receiver  are  at 
least  several  wavelengths  from  the  sea  surface,*  ray  theory  can  be  accurate  in  its  predictions  for  these 
ducted  environments. 

Study  Motivation 


This  study  is  concerned  with  assessing  the  capabilities  and  limitations  of  ray-theoretic  models 
when  used  to  predict  mid-  and  long-range  signal  behavior  in  deep-water  environments  with  superim¬ 
posed  surface  ducts.  The  ray-theoretic  model  is  one  which  has  been  incorporated  as  a  submodel  into 
the  NRL  reverberation  model  [2],  The  PE  model  used  for  this  study  is  one  that  has  also  been 
developed  at  NRL  [3], 

In  general,  the  PE  method  (Appendix  A)  has  achieved  widespread  acceptance  in  the  acoustics 
community  because  of  its  versatility  over  a  large  frequency  domain  and  under  many  environmental 
conditions.  However,  such  an  approach  is  not  always  appropriate  or  practical,  particularly  when  the 
time-dependence  of  a  pulse  signal  and  when  backscattered  energy  must  be  considered  as  in  reverbera¬ 
tion  studies.  An  alternate  method  often  used  is  the  ray-theoretic  approach  (Appendix  B).  Ray  theory 
has  a  number  of  advantages  over  the  PE  method  for  computing  reverberation  in  long-range  active  sonar 
applications.  These  advantages  are: 

•  It  easily  computes  time  of  arrival  for  pulsed  signals. 

•  Models  for  backscattering  and  forward  scattering  of  energy  can  be  easily  incorporated. 


There  is  no  limitation  on  the  propagation  angles  which  can  be  considered. 


Manuscript  approved  March  31,  1982. 

•When  source  and  receiver  are  at  the  surface,  an  unmanageable  piling  up  of  caustics  as  range  increases  results  in  the  complete 
breakdown  of  ray  theory  [1],  The  results  show  no  such  difficulties  for  a  source  at  91  m  receiver  of  20  m  to  ranges  of  150  nmi. 


1 


TOLSTOY,  FRANCHI,  AND  NICOLAS 


The  propagation  model  used  in  this  study  is  a  corrected  ray-theoretic  model  which  overcomes  the 
normal  difficulties  with  caustics  and  discontinuities  in  derivatives  of  the  sound-speed  profile.  The 
model  incorporates  a  correction  for  smooth  caustics  [4]  based  upon  a  uniform  asymptotic  expansion  for 
the  inverse  Hankel  transform  of  the  WKB  solution  of  the  Hankel  transformed  wave  equation.  The 
solution  can  be  extended  into  the  caustic  shadow  zone  to  predict  appropriate  nonzero  intensities.  More¬ 
over,  the  sound-speed  profile  data  points  are  fit  by  a  cubic  spline  interpolating  polynomial  to  ensure 
good  behavior  for  the  first  and  second  depth  derivatives. 

The  ray-theoretic  model  can  be  expected  to  have  limitations  at  low  frequencies,  under  conditions 
of  strong  sound-speed  gradients,  and  when  diffraction  and  mode  coupling  effects  are  dominant.  How¬ 
ever,  the  extent  to  which  these  disadvantages  limit  the  applicability  of  a  medium-  to  long-range  rever¬ 
beration  model  incorporating  a  corrected  ray-theoretic  propagation  submodel  needs  to  be  addressed. 
This  study  was  motivated  by  the  belief  that  a  ray-based  reverberation  model  could  provide  reasonable 
calculation  results  for  long-range  active  sonars  under  a  variety  of  important  environmental  conditions, 
including  strong  winter  surface  ducts.  This  study  is  aimed  at  assessing  such  a  capability  by  comparing 
the  propagation  results  from  the  two  models  for  environmental  conditions  in  which  the  PE  method  had 
been  used  in  a  previous  study  and  for  which  those  results  have  received  wide  acceptance.  Our  study 
uses  the  same  sound-speed  environments  as  the  previous  PE  studies,  but  only  some  of  the  source- 
receiver  placements.  Thus,  this  study  is  not  an  exhaustive  investigation,  but  one  in  which  a  number  of 
conditions  of  general  interest  are  explored.  The  objective  is  to  see  how  well  a  ray-theoretic  propagation 
model  will  work  in  situations  where  it  would  be  practical  to  use  such  a  model  as  part  of  a  reverberation 
prediction  capability.  Frequently  such  ray-based  approaches  to  reverberation  estimation  are  dismissed 
without  quantitative  justification  even  though,  they  can  achieve  good  accuracy  and  save  computer  time. 

TEST  DATA 

Three  ducted  environments  selected  from  measured  data  have  been  considered  for  this  study. 
Their  surface  ducts  range  from  200  m  to  800  m  in  depth  while  the  total  water  depth  is  a  constant  3000 
m.  Although  the  ray  program  can  handle  variations  of  bottom  topography  and  sound  speed  with  range, 
we  elected  to  examine  only  flat  bottom  and  range-independent  environments.  In  this  manner  we  could 
most  confidently  deduce  the  causes  of  any  propagation  discrepancies.  See  Fig.  1  and  Table  1  for  the 
profiles  selected  for  this  study  and  the  comparison  of  properties  of  the  sound-speed  profiles.  The 
source  depths  selected  were  91  m  and  305  m  while  receiver  depths  were  at  20,  50,  90,  and  306  m;  we 
also  presume  a  point  source  with  a  300-Hz  constant  frequency.  For  waveguide  boundaries  we  assume  a 
totally  absorbing,  liquid  bottom  and  a  flat,  totally  reflecting  (with  phase  change  of  tt)  sea-surface.  Note 
that  the  shallowest  duct  considered  here  is  still  some  40  wavelengths  deep,  and  the  maximum  gradient 
dc/ dz  computed  for  all  profiles  was  about  0.027  s'"1.  Thus,  0.027/300  «  4  tt  implies  that  ray-theory 
applicability  is  not  contraindicated. 

RESULTS 

Figures  2a  to  7a  represent  PE  signal  intensity  plots  for  profiles  1,  2,  and  3  as  indicated  with  source 
depths  located  at  91  m  and  305  m.  Figures  2b  to  7b  represent  the  corresponding  ray  traces.  In  the  PE 
figures,  the  pronounced  and  genuine  (measurable)  interference  patterns  (the  shading  with  dark  regions 
indicating  high  intensity).  These  patterns  are  produced  by  the  coherent  addition  of  the  PE  modes/rays. 
The  ray  figures  display  more  diagrammatic  patterns  or  ray  tracks  since  the  intensity  calculations  are  per¬ 
formed  at  a  later  stage  in  the  program’s  computations.  However,  a  number  of  conclusions  may  still  be 
drawn  by  comparison  between  these  ray  and  PE  plots.  First,  for  profiles  1  and  2  (Figs.  2  to  5)  both 
types  of  plots  indicate  significant  levelS/Of  energy  trapped  in  the  surface  duct.  Next,  they  both  indicate 
the  presence  of  deep  cycling  energy  which  has  escaped  the  surface  duct  and  which  displays  decreasing 
focusing  (more  smearing  out)  with  increasing  range  from  the  sources.  Examining  the  behavior  elicited 
by  profile  3  (Figs.  6  and  7)  we  note  that  for  this  shallower  surface  duct  (200  m  deep),  the  source  at  305 
m  is  located  below  the  duct.  As  a  result,  little  signal  energy  becomes  trapped  near  the  surface.  Finally, 
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Fig.  1  —  Three  ducted  sound-speed  profiles.  Data  points  are  indicated 
by  points  with  spline  interpolation  used  between  them. 


Table  1  —  Comparison  of  Properties  of  Sound-Speed  Profiles 


Property 

pi 

P2 

P3 

Surface  duct  depth  (m) 

800 

381 

200 

SOFAR  axis  depth  (m) 

1334 

1334 

724 

Ocean  depth  (m) 

3000 

3000 

3000 

Sound  speed  at  the  ocean  surface  (m/s) 

1483 

1496 

1487 

Sound  speed  at  surface  duct  maximum  (m/s) 

1495 

1501 

1490 

Sound  speed  at  SOFAR  axis  (m/s) 

1490 

1490 

1482 

Sound  speed  at  ocean  bottom  (m/s) 

1512 

1512 

1512 

we  see  that  all  the  PE  figures  agree  extremely  well  with  the  ray  figures  in  the  fine  detail  of  the  signal 
structure,  and  indeed,  if  these  plots  are  overlayed,  it  becomes  apparent  that  the  patterns  are  in  excellent 
alignment  with  each  other.  To  understand  conditions  under  which  noticeable  shifts  in  such  patterns 
might  occur  (between  the  ray  and  PE  approaches)  we  also  considered  the  effects  of  two  strongly  refract¬ 
ing  analytic  sound-speed  profiles.  Appendix  C  contains  these  results. 
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Fig.  2  —  Ray  tracks,  PE  intensity  plots  and  transmission  loss  plots  as  a  function  of  range 
for  sound  speed  profile  1  with  SD  —  91  m,  RD  *  305  m,  300  Hz 
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Fig.  3  —  Ray  tracks,  PE  intensity  plots  and  transmission  loss  plots  as  a  function  of  range 
for  sound  speed  profile  1  with  SD  =  305  m,  RD  —  20  m,  300  Hz 
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Fig.  4  —  Ray  tracks,  PE  intensity  plots  and  transmission  loss  plots  as  a  function  of  range 
for  sound  speed  profile  2  with  SD  =  91  m,  RD  =  20  m,  300  Hz 
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Fig.  5  —  Ray  tracks,  PE  intensity  plots  and  transmission  loss  plots  as  a  function  of  range 
for  sound  speed  profile  2  with  SD  =  305'm,  RD  =**  20  m,  300  Hz 
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Fig.  6  —  Ray  tracks,  PE  intensity  plots  and  transmission  loss  plots  as  a  function  of  range 
for  sound  speed  profile  3  with  SD  =  91  m,  RD  =  50  m,  300  Hz 
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Fig.  7  —  Ray  tracks,  PE  intensity  plots  and  transmission  loss  plots  as  a  function  of  range 
for  sound  speed  profile  3  with  SD  —  305  m,  RD  —  20  m,  300  Hz 
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Next,  we  examined  predicted  levels  of  signal  transmission  loss  vs  range  (at  a  specified  receiver 
depth).  These  levels  are  plotted  in  both  raw  and  range  averaged  (Gaussian  weights  over  6  cr  —  5  km) 
forms.  In  Figs.  2c-7c  and  8a  the  ray  model  plots  now  represent  both  an  incoherent*  (dashed  line)  and  a 
coherent  sum  of  ray  intensities.  We  can  expect  some  disagreement  between  the  two  models  at  ranges 
near  the  source  (within  20  km)  because  the  PE  approximation  is  expected  to  acquire  validity  only  at 
longer  ranges  (we  recall  the  far-held  and  paraxial  assumptions).  In  addition,  for  the  ray  model  the  den¬ 
sity  of  rays  has  been  limited  deliberately  at  the  higher  (more  vertical)  launch  angles  which  can  some¬ 
times  result  in  unrealistic  zero  intensities  (large  TL)  in  the  short  range  unaveraged/raw  data  (see  Fig.  8 
in  the  ranges  5  to  10  km).  This  sparsity  should  not  be  regarded  as  more  than  a  temporary  and  inconse¬ 
quential  aberration.  That  is,  the  high-angle  rays  which  would  be  included  for  a  more  accurate  near-held 
effect  eventually  interact  often  enough  with  the  bottom  to  become  absorbed  and  thereby  not  affect 
results  at  the  longer  ranges  which  are  the  focus  of  this  study.  For  prohle  3,  source  depth  305  m, 
receiver  depth  20  m  we  again  note  zero  intensity  levels  in  the  ray  model  predictions  (Fig.  7d).  How¬ 
ever,  we  note  here  that  the  gap  is  rather  extended  in  range  (<  20  km)  and  reappears  several  times 
(ranges  13  to  38,  60  to  80,  113  to  125  km)  before  the  deep  cycling  energy  smears  out  sufficiently  near 
the  surface  to  introduce  consistently  nonzero  intensity  levels.  These  gaps  (for  a  receiver  in  the  duct, 
source  below  the  duct  or  vice-versa)  are  due  primarily  to  the  failure  of  the  ray  theoretic  approach  to 
consider  diffracted  energy.  Note,  however,  that  these  gaps  would  have  been  partially  filled  if  bottom- 
reflected  energy  had  been  considered  (a  more  realistic  situation). 

We  conclude  that  the  levels  of  transmission  loss  vs  range  predicted  by  the  two  models  are  in  good 
agreement  (see  Figs.  2c,  2d,  7c,  7d,  8).  Both  models  involve  different  theoretical  as  well  as  numerical 
approximation  to  the  exact  solution.  Consequently,  we  should  not  expect  the  detailed  interference  pat¬ 
terns;  that  is,  the  location  of  signal  maxima  and  minima,  as  predicted  by  each  model  to  agree,  particu¬ 
larly  as  range  increases.  In  the  raw  data,  where  the  oscillations  are  quite  frequent,  a  statistical  compari¬ 
son  would  be  required  to  show  how  well  the  two  models  agree  or  disagree  with  regard  to  fine  structure. 
It  is  to  be  expected  that  the  PE  and  ray-theoretic  approaches  should  not,  in  fact,  yield  the  exactly  ident¬ 
ical  responses  to  any  given  environment.  Indications  of  when  and  how  the  two  solutions  can  differ  may 
be  found  in  Appendix  C  where  two  analytical  profiles  are  considered. 

CONCLUSIONS 

For  the  test  cases  examined,  the  ray  model  and  the  PE  model  agree  extremely  well.  In  particular, 
the  ray  tracks  and  PE  intensity  plots  show  excellent  and  detailed  correlation  in  the  ducted-  and  deep¬ 
cycling  energy.  The  range-averaged  TL  levels  also  agree  except  under  some  near-field  conditions  where 
either  the  lack  of  a  sufficient  number  of  high-angle  rays  or  the  lack  of  consideration  of  diffraction 
effects  in  the  ray  model  has  occasionally  led  to  levels  which  are  unrealistically  low.  Thus,  the  few  occa¬ 
sions  on  which  the  two  approaches  differed  are  well  understood.  Since  the  focus  of  this  study  is  on 
long-range  prediction,  these  areas  of  disagreement  are  not  considered  serious.  As  a  result,  we  believe 
that  long-range  predictions  as  made  by  the  extended  ray  model  used  in  this  report  should  be  as  accurate 
as  PE  predictions  in  ducted  environments  with  \  <  h/\0  and  for  sources  and  receivers  away  from  the 
sea-surface. 


This  is  essentially  a  weighted  (Gaussian  in  range  and  centered  about  the  receiver  location)  sum  of  ray  intensities  at  each  depth 
12,51. 
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Fig.  8  —  Unaveraged  transmission  loss  as  a  function  of  range  (a:  ray  theoretic;  b:  PE  theoretic) 
for  sound  speed  profile  1  with  SD  =  91  m,  RD  =  305  m,  freq  =  300  Hz. 
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Appendix  A 
PE  BACKGROUND 


By  linearizing  both  the  Navier-Stokes’  equation  for  a  nonviscous  fluid  (Newton’s  equation  of 
motion),  and  the  equation  of  continuity  (conservation  of  mass),  one  finds  that  to  a  good  first  approxi¬ 
mation,  propagation  of  small  amplitude  sound  waves  in  an  ocean  environment  can  be  described  by  the 
linear,  hyperbolic,  second-order,  time-dependent,  scalar  wave  equation  [A1,A2]: 

V20 - y~—  =  -f(x,t),  (Al) 

czOc)  or 

where  ifj  is  the  acoustic  pressure  perturbation, 

x  is  the  position  vector  in  a  3-dimensional  coordinate  system, 
c(x)  is  the  sound  speed  profile  (assumed  time-independent),  and 
f(x,t)  is  the  source  function. 

By  using  appropriate  boundary  and  initial  conditions  (BCs  and  ICs)  one  is  guaranteed  the  existence  of  a 
unique  solution  in  some  appropriate  solution  space.  However,  as  c(x)  or  the  BCs  become  complicated, 
it  becomes  impossible  to  compute  an  exact  analytic  solution.  Hence,  schemes  to  approximate  this  solu¬ 
tion  have  been  developed,  one  of  which  is  the  parabolic  approximation. 


Let  us  now  consider  a  two-dimensional  environment  with  cylindrical  coordinates  ( r,z )  where  r  is 
range,  z  is  depth  (z  =  0  at  sea  surface,  z  >  0  below  surface).  Then  for  harmonic  time  dependence  the 
wave  equation  becomes  the  Helmholtz  (reduced)  wave  equation: 


V2t jj(r,z)  +  k$n2(r,z)i(j(rtz)  =  —f(r,z) 


(A2) 


where  /c0  is  the  reference  wave  number, 

rf(i)z)  is  jthe  index  of  refraction, 

■\«/ ( k0(x )),  and 
a)  is  the  source  angular  frequency. 


We  can  write  the  Green’s  function  solution  (/  =  8  (r)  6  (z  -  z0))  as 


iknr 


0  ( r,z )  =  — j=~  <f>(r,z ), 


where  </>  satisfies  (away  from  the  source): 


HA  +  o ik  +  HA  +  l.2 

2  +  UKq  +  +  k0 

or  or  oz 


If  we  now  assume  that 


and  that 


92<ft 


dr2 


« 


n2(r,z)  -  1  + 


d<t> 


4kb2 


<f>  =  0. 


2k°  dr 


=  0  (far  field), 


rr 


then 


2ik0  +  HA  +  k 2  ( n2(r,z )  —  1)0  ~  0. 


dr  dz 2 


(A3) 

(A4) 

(A5) 

(A6) 

(A7) 
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Equation  (A5)  is  usually  referred  to  as  the  paraxial  approximation  which  implies  that  if  the  acoustic 
field  were  represented  by  rays  then  these  would  be  included  only  at  small  angles  with  respect  to  the 
horizontal  [A3-A5]. 


A  solution  to  Eq.  (A7)  may  be  found  simply  by  Fourier  transforming  the  equation  (on  z)  and 
then  assuming  that 


then  integrating  by  parts  to  obtain 


where 


Then, 


n{r,z)  =  constant, 

<f>,  — 1 *■  0  as  z  -+ 

oz 


3$  ,  k§  ( n 2  -  1)  -  s2  ^ 

-7—  + - 777 - *  =  0, 

or  2  ik0 


def  i  /» oo 

<t >(r,s)  —  - —  I  <j>(r,z)eiszdz. 
2  77  J-°° 


94> 

dr  s2  —  ko  ( n 2  —  1) 


$ 


2ikn 


So,  integrating  from  rto  r  +  A 7  we  have 

In 

implying 


<J>(r  +  A  r,s) 

•  =  A  r  • 

V- 

ko  ( n 2  -  1) 

3>(r,s) 

2  ik0 

d>(r  +  A r,s)  —  <&(/•, s)  exp 


This  equation  can  be  inverse  transformed  to  obtain 


iAr 

2 


k0(rt2  -  1)  — j— 
ko 


,  /  ,  .  \  ikn(n2  -  1)  A T  ru-i 

0(r  +  Ar,z)  =  e  0  — 


A  r 


2k, 


D(r,s) 


where 


def 


J^M/}  =  f"  f{s)e~iszds. 

4/  — oo 


(A8) 

(A9) 

(A10) 

(All) 

(A12) 

(A13) 

(A  14) 

(A  15) 

(A  16) 


Thus,  the  field  at  any  range  ( r  +  Ar)  depth  zis  obtained  by  stepping  the  solution  by  a  range  increment 
Arfrom  the  earlier  value  at  range  r,  depth  z. 


Theoretically  and  numerically  this  solution  is  stable  and  highly  accurate  even  when  n(r,z)  is  not 
constant  provided  Ar  is  sufficiently  small  [A4],  Within  the  conditions  of  its  application,  this  treatment 
includes  all  effects  normally  associated  with  diffraction  and  mode  coupling  in  a  range  and  depth  depen¬ 
dent  environment,  but  does  not  allow  for  the  backscattering  of  acoustic  energy  or  for  the  propagation  of 
high  angle  rays.  The  sea-surface  (z  =  0)  boundary  condition  of  perfect  reflection  plus  a  phase  change 
of  77-  (0  =  0  at  z  =  0)  is  replaced  by  an  image  source  at  z  =  — zO.  Also,  for  |z|  >  bottom  depth,  n2  is 
given  an  exponentially  increasing  imaginary  part  resulting  in  a  solution /which  damps  to  zero  as 
|z|  — • ►  oo.  An  initial  field  ( r  =  0)  must  be  computed  before  the  marching  algorithm  (split-step)  can 
begin  operating  on  Eq.  (A  15),  and  this  field  has  been  calculated  here  as  a  very  narrow  Gaussian  field 
for  the  assumed  point  source. 
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The  PE  model  used  [A5]  has  a  number  of  attractive  features. 

•  All  frequencies  can  be  accurately  modeled  by  suitably  choosing  A  r. 

•  A  range  variable  (as  well  as  depth  variable)  environment  can  be  considered. 

•  Diffraction  and  mode  coupling  effects  are  implicitly  included. 

•  A  choice  of  initial  field  is  allowed,  e.g.,  Gaussian,  normal  mode,  user  specified. 

•  Solution  is  highly  stable. 

•  Sound  speed  profiles  need  only  simple  linear  interpolation  in  depth  between  data  points. 
However,  a  PE  approach  also  has  its  limitations. 

•  Time  of  arrival  for  pulsed  signals  (finite  duration)  cannot  be  computed. 

•  Backscattered  energy  is  not  allowed. 

•  The  paraxial,  that  is  small  angle,  approximation  must  be  assumed  (an  acceptable  assump¬ 
tion  for  most  long  range  studies). 

•  Program  requires  large  amounts  of  file  storage  for  long  range/ high  frequency  calculations. 
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Appendix  B 

RAY  THEORY  BACKGROUND 


The  geometrical  approximation  can  be  derived  in  a  number  of  ways.  One  approach  (we  shall  hen¬ 
ceforth  assume  a  harmonic  time  dependence)  is  to  assume  [Bl]  that  the  solution  must  be  of  a  certain 
form,  that  is, 


<K3c)  =  A(x)ei(k°S™-ml) 


(Bl) 


where 


a)  is  the  source  angular  frequency, 
k0  is  the  reference  wave  number, 
AC x)  is  the  amplitude  at  x,  and 
Six)  is  the  phase  at  x. 


Substituting  this  expression  into  the  homogeneous  wave  equation  fix,t )  =  0  in  Eq.  (Bl)  and  col¬ 
lecting  real  and  imaginary  terms  leads  to 

V2A  -  fc02  iS2  +  S2  +  S2)  A  +  -j- ]-  A  =  0,  and  (B2) 

C  (  X) 

2 V A  •  VS  +  A  V2S  =  0  (transport  equation).  (B3) 

If,  we  now  assume  that 

<K  =  (B4) 

where  the  index  of  refraction  n  (x)  is 


nix)  = 


OJ 


k0c(x )  ’ 

then  Eq.  (B2)  becomes  the  standard  eikonal  equation  (independent  of  source  frequency): 

|VS|2=  n2. 

We  note  that  inequality  Eq.  (B4)  may  not  hold. 

•  At  too  low  frequencies  («  small),  or 

•  for  V2A  large. 


(B5) 


The  computation  of  signal  amplitude  A  (and  subsequently,  of  V2A)  in  ray  theory  proceeds  either  by 
solving  Eq.  (B3)  when  possible  or  by  invoking  the  principle  of  conservation  of  energy  flux  for  a  ray 
bundle.  If  we  assume  constant  density  this  leads  to  [Bl] 


A  (x2) 


ko 


cix2) 

ciXi) 


dai 
d(T  2 


1/2 

A 


ix\), 


(B6) 


where 


A  (%)  is  the  field  amplitude  at  xh  i  =  1,  2 
d<T\  is  the  ray  bundle  cross  section  at  x,,  i  =  1,2  (Fig.  Bl). 
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Fig.  B.l  —  Calculation  of  length  dLi  of 
ray  bundle  cross  section  at  xf 


Clearly,  Eq.  (B6)  predicts  A  (x2)  —  (for  A  (xi)  <  oo,  dv\  <  oo)  when  da2  0-  That  is,  ray 
theory  does  not  hold: 

near  caustics  and  focal  points. 


On  the  basis  of  geometric  arguments  (Fig.  Bl)  we  have  for  a  point  source  that 

do-j  =  IttRi  sin0,  dRif  (Bl) 

where  Rf  is  the  horizontal  distance  from  source  to  xh  and 
0,  is  the  ray  angle  (w.r.t.  horizontal)  at  xh 

Now,  when  0  =  0,  that  is  at  the  ray  turning  point,  sin  0  — ■ ►  0  and  dR  In  this  instance,  however, 

the  product  has  a  well-defined  limit,  that  is,  sin  QdR  =  cos  Odz  =  dz  when  0  — ►  0.  Thus,  for  a  point 
source  ray  theory  gives  accurate  amplitude  estimates  at  turning  points.  This  is  not  true  for  a  "plane- 
wave"  (WKB)  solution  [Bl]  where  the  turning  points  lie  on  a  caustic. 

To  obtain  an  alternate  expression  for  A  (x2)  let  xx  be  at  unit  radial  distance  from  source,  that  is, 
at  R\  —  cos  0i  where  0X  is  the  ray  launch  angle.  Let  cGq)  be  the  reference  sound  speed,  Aq  =  A  (xi) 
the  reference  (source)  amplitude.  Then,  dzx  cos  0!  =  dO  implies  that 


A  (x2) 


rn 

V  n(x2) 


cos d6 


«(xi)  Ri  sin  0 2  dRi 


An. 


One  final  cautionary  note  on  the  limits  of  the  validity  of  ray  theory  must  also  be  presented,  and 
this  is  most  easily  developed  within  the  context  of  a  depth  stratified  medium  (cCx)  =  c(z)). 


We  notice  that  A  (x)  —  V c(x)  =  Vc(z);  so,  if  we  now  also  assume  that  we  are  away  from  turn¬ 
ing  points  then  V2A  ~  Azz  and 


Azz  ~  -  {c(z)}1/2  •  A 

dz 


(B8) 


turning  point 


1 

c- 1/2  _ 

,  / 

_L  r-3/2 

dcX 

r 

•  A 

2 

dz1 

2 

dz 

turning  point 
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(p-  Q 

Thus  ~7J  sma11  ,eads  t0  (with  A  =  VcIzT  A  |turnin„  D0int) 


dz 


-  1  -2 

dc 

A 

4  C 

dz 

Therefore, 


V2 * 

1  r~2 

dc 

2 

(JL>  2 

A 

4  C 

dz 

« 

c2 

i 

dc 

2oj 

dz 

«  1. 


(B9) 


(BIO) 


This  inequality  (BIO)  illustrates  an  intuitive  low-gradient  requirement  (necessary  but  not  sufficient)  for 
the  ray  approximation  [B1,B2], 


A  wave-theoretical  approach  leading  to  an  optical  theoretic  solution  is  to  construct  the  Green’s 
function  in  integral  form,  reexpress  the  function  as  a  series  in  terms  of  traveling  waves,  and  then 
reduce  the  integrals  by  means  of  saddle-point  evaluations  [B3,B4], 


The  present  ray-theoretic  model  expands  the  parameters  of  interest  for  each  ray,  that  is,  range, 
depth,  angle,  travel  time,  Snell’s  "constant,"  by  means  of  a  Taylor  series  expansion  (finite  number  of 
terms)  expressed  as  a  function  of  arc-length,  /.  These  terms  are  evaluated  by  means  of  the  ray  equa¬ 
tion: 


d_ 

dl 


dx 

dS 


=  Vn, 


(Bll) 


where  /  is  the  arc  length  along  the  ray,  and  by  means  of  environmental  data,  that  is,  first-  and  second- 
order  derivatives  of  c  ( x )  with  respect  to  depth  and  range. 

Equation  (Bll)  can  be  derived  rigorously  from  Fermat’s  principle  (the  ray  is  path  of  minimum 
travel  time  between  two  fixed  points)  by  the  calculus  of  variations.  It  can  also  be  derived*  from  the 
eikonal  equation  by  noting  that  the  unit  vector  normal  p  to  the  wavefront  (surface  of  constant  phase)  is 
given  by 


P  = 


VS 


(follows  from  VS 


=  n2)  and  is  also  given  by  the  unit  ray  tangent  v  where 


(B12) 


„  dx  i  dx2  dx  3 

V=  [dT’  ~df’  ~di 

where  (xj(/),  x2(l),  x2(l))  is  point  on  the  ray.  Thus, 

dx,  l  QS 

dl  n  dx/  ’ 


(B13) 


(14) 


*This  derivation  is  courtesy  of  D.  Berman  (NRL). 
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Differentiating  by  /  we  get 


d_ 

dl 


dx, 

~dl 


d_  dS_ 
dl  dXj 


d 

dx, 

d 

dx, 

_d_  l 

dx, 

d 

dx, 

d 

dx, 


dS 


dl 

dS  dx i 


dx 

ivs  •  V 


_  _dS_  dxi_  ,  dS_  <&i 

dl  dx2  dl  9x3  dl 


VS-p\ 

( n ) 


since 


V  S  =  np. 


(B15) 


The  model  used  in  this  study  assumes  azimuthal  symmetry  for  the  environment  and  hence,  it 
represents  an  approximate  solution  to  the  two-dimensional  rather  than  to  the  full  three-dimensional 
wave  Eq.  (Al)  in  Appendix  A  of  this  report. 

In  general  ray  theory  has  a  number  of  advantages. 

•  It  easily  computes  time  of  arrival  for  pulsed  signals. 

•  It  is  versatile  enough  to  handle  complex  geophysical  conditions. 

•  Models  for  backscattering  and  forwardscattering  of  energy  can  be  easily  incorporated. 

•  It  is  an  easily  visualized  technique. 

•  The  equations  are  more  tractable  than,  say,  normal-mode  equations,  in  a  realistic  environ¬ 
ment. 

•  There  is  no  limitation  on  the  propagation  angles  which  can  be  considered. 

Of  course  ray  theory  also  has  its  disadvantages. 

•  It  does  not  predict  well  at  low  frequencies,  for  example, 
when  X  > 


where  X  is  signal  wavelength,  and 
h  is  duct  depth  [B5], 
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or  under  strong  gradient  conditions  when 


_1_ 

2ct> 


dc 

dz 


>  1, 


where  a>  is  angular  source  frequency, 

c(r,z)  is  the  sound  speed  profile, 
r  is  horizontal  range, 
z  is  depth. 


The  determination  of  appropriate  rays  (those  connecting  source  and  receiver)  can  be 
difficult  or  inefficient. 


•  The  theory  needs  modification  near  caustics  and  focal  points. 

•  Its  predictions  are  very  sensitive  to  discontinuities  in  the  slopes  and  curvatures  of  the 
sound  speed  profiles  [B6-B8], 

The  model  used  an  automated  ray-selection  procedure  which  can  be  supplemented  or  replaced  by 
user  denned  rays. 
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Appendix  C 

ANALYTIC  PROFILE  RESULTS 

It  can  be  shown  [C1,C2]  that  a  ray-theoretic  solution  should  focus  periodically  and  independently 
of  launch  angle  when  the  sound-speed  profile  is  of  the  form 

c(z)  =  A  cosh  Biz  —  C), 

and  when  the  source  and  receiver  are  located  at  the  profile  axis,  that  is,  at  z  =  C.  For  this  test  case  we 
selected  A  =  1500  m/s,  B  =  0.003  m-1,  and  C  =  1500  m  for  a  total  water  depth  of  3000  m  (Fig.  Cl). 
Figure  C2a  shows  how  well  the  ray  model  reproduces  the  predicted  focusing.  In  contrast,  Fig.  C2b 
shows  the  lack  of  perfect  focusing  by  the  PE  model  for  this  profile.  Next,  we  considered  a  profile  for 
which  the  PE  is  predicted  to  focus  perfectly,  that  is  one  for  which  the  square  of  the  index  of  refraction 
is  parabolic.  We  matched  this  profile  as  closely  as  possible  to  the  previous  hyperbolic  cosine  profile, 

c(z)  =  =  — 

a/— /3(z  —  1500)2  +  a  ’ 

where  c(0)  =  c(3000)  =  1654.4  — »  a  =  1 

j8  =  7.9x  10"8 

(V/3  =  2.8  x  10-4). 

SOUND  SPEED  (M/S) 
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Note  that  if 
fix) 


and 

g(x) 


(a)  Ray 


RANGE  (KM) 
(b)  PE 


PROFA  (COSH)  PE 


(DB/ZM) 


Fig.  C.2  -  Ray  tracks  and  PE  intensity  plots  for  the  hyperbolic 
cosine  sound-speed  profile  with  freq  =  300  Hz 


cosh  bx 

ebx  +  e-bx 

2 

~  Y  {2  +  ( bx )2} 


y  ( bxYn  +  (-to)" 
„Tn  n\ 


1  +  (1/2)  (to)2. 


1 

VI  -  /to2  ’ 
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then  . 

fix)  =  gix )  for  b2  ^  ft. 

We  have  for  the  above  examples  V/3  =  2.8  x  10~4,  fc  =  3.0x  10-4  (Fig.  Cl).  Then,  in  Fig.  C3b  we 
see  the  PE  focusing  perfectly  while  in  Fig.  C3a  the  ray  model  does  not.  This  shift  is  also  apparent  in 
the  TL  vs  range  plots  at  a  receiver  depth  of  100  m  (Fig.  C4a,b).  Hence,  we  conclude  that  each  model 
can  react  differently  to  the  same  sound-speed  environment  and  show  offsets  between  their  respective 
interference  patterns.  Moreover,  there  is  very  little  difference  between  these  two  analytical  profiles,  but 
that  difference  is  sufficient  to  produce  noticeable  disagreements  in  the  model  predictions;  yet  the 
models  are  each  performing  exactly  as  they  should.  As  a  result,  we  should  not  expect  perfectly  corre¬ 
lated  TL  plots  between  the  models  for  any  given  environments.  Note  that  the  analytic  profiles  selected 
here  were  much  more  refracting  than  one  is  likely  to  see  in  a  real  environment  (min  c(z )  =  1500.0 
m/s,  max  ciz)  =  1654.5  m/s),  and  as  such  the  model  differences  have  been  much  more  evident  for 
these  cases  than  one  is  likely  to  ever  encounter  in  a  realistic  test. 


TEST  PROFILE  (PARABOLIC  N»m2)  -Ml*-  Jd*  U t 


(a)  Ray 

PROFB  ( N**2  PRRfiBOL I C )  PE 


O  40 

RfiNGE  (KM) 

(b)  PE 


Fig.  C.3  —  Ray  tracks  and  PE  intensity  plots  for  the  n2  parabolic 
sound-speed  profile  with  freq  =  300  Hz 
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TEST  PROFB 


n2U)  PARABOLIC 
SD= 1500  m 
RD=  100  m 
FREO 1 300  Hz 
— Rot— 

(a) 


n2(z)  PARABOLIC 
SD  B 1 500  m 
RD=  100  m 
FREQ « 300  Hz 


Fig.  C.4  —  Unaveraged  transmission  loss  as  a  function  of  range  (a:  ray  theoretic,  b:  PE  theoretic)  for 
the  n 2  parabolic  sound-speed  profile  with  SD  «  1500  m,  RD  =  100  m,  freq  =  300  Hz 
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